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Abstract 

We consider nonleptonic B decays into two light mesons at leading order in soft- 
collinear effective theory, and show that the decay amplitudes are factorized to all 
orders in a s . The operators for nonleptonic B decays in the full theory are first 
matched to the operators in SCETi, which is the effective theory appropriate for 
\Art&A < fjL < mi, with A ~ 0.5 GeV. We evolve the operators and the relevant 
time-ordered products in SCETi to SCETn, which is appropriate for u. < ^/m^K. 
Using the gauge-invariant operators in SCETn, we compute nonleptonic B decays 
in SCET, including the nonfactorizable spectator contributions and spectator con- 
tributions to the heavy-to-light form factor. As an application, we present the decay 
amplitudes for B — ► tttt in soft-collinear effective theory. 



1 Introduction 



Decay of B mesons plays an important role in particle physics since it is a test- 
ing ground for the Standard Model and a window for possible new physics. 
We can obtain good information on the Cabibbo-Kobayashi-Maskawa (CKM) 
matrix elements and CP violation from B decays. Since the contribution of 
QCD in B decays changes the whole structure of the theory, the study of B 
decays is an intertwined field of particle physics. Among these decays, non- 
leptonic B decays have been the subject of intense interest. Especially the 
treatment of nonperturbative effects from the strong interaction is a serious 
theoretical problem in nonleptonic decays. Precise experimental observation 
of nonleptonic B decays makes it urgent to give firm theoretical prediction in- 
cluding the effects of CP violation. There has been a lot of theoretical progress 
in nonleptonic B decays and we suggest how to consider nonleptonic B decays 
from the viewpoint of the soft-collinear effective theory (SCET). 
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The effective Hamiltonian for nonleptonic decays from the Standard Model 
has been derived and the Wilson coefficients of the operators for B decays 
have been calculated to next-to-leading order, and next-to-next-to-leading or- 
der for some operators [1] . In order to evaluate the hadronic matrix elements of 
four-quark operators, naive factorization [2] was assumed, in which the matrix 
elements were reduced to products of current matrix elements. But there was 
no justification for this assumption except the argument of color transparency 
[3]. Besides that, decay amplitudes depend on an arbitrary renormalization 
scale pi in naive factorization since the Wilson coefficients depend on //, while 
the matrix elements of operators do not. Ali et al. [4] have improved the prob- 
lem of the scale dependence by including radiative corrections of the operators 
before taking hadronic matrix elements. Then the \x dependence of the Wilson 
coefficients is cancelled by that of the radiative corrections of the operators. 
However, the decay amplitudes depend on calculation schemes since the off- 
shell renormalization scheme is used [5]. 

Politzer and Wise [6] suggested to take the heavy quark mass limit in comput- 
ing corrections to the decay rate ratio T(B — > D*rr)/T(B — > Dn). They have 
considered radiative corrections for nonfactorizable contributions and found 
that they are finite and the decay amplitudes are factorized. Beneke et al. [7] 
have extended this idea to general two-body decays including two light final- 
state mesons, in which the decay amplitudes can be expanded in a power series 
of 1/mf,. They show that nonfactorizable contributions including spectator in- 
teractions are factorized as a convolution of the hard scattering amplitudes 
with the meson wave functions, and the corrections are suppressed by pow- 
ers of 1/rrib. This is an important step toward theoretical understanding of 
nonleptonic B decays. First the amplitudes can be obtained from first princi- 
ples and a systematic 1/m^ expansion is possible. Second, since the on-shell 
renormalization is used, there is no scheme dependence. And they improved 
previous approaches by including momentum- dependent parts, which had not 
been included previously. However, when higher-twist light-cone wave func- 
tions are included, there appear infrared divergences in the amplitudes. These 
are treated as theoretical uncertainties, but from the theoretical viewpoint it 
is a problem to be solved in this approach. 

Bauer et al. [8,9] have proposed an effective field theory in which massless 
quarks move with large energy. This effective theory, called the soft-collinear 
effective theory, is appropriate for light quarks with large energy. It has been 
applied to hard scattering processes and B decays [10-15]. It is also an ap- 
propriate effective theory for nonleptonic B decays to two light mesons. In 
this paper, we apply SCET to nonleptonic B decays into light mesons. We 
construct all the relevant operators in the effective theory at leading order 
in a gauge-invariant way by integrating out all the off-shell modes. The Wil- 
son coefficients are calculated by matching the full theory onto SCET. We 
show that the four-quark operators in SCET are factorized to all orders in a s 
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and the argument of color transparency is explictly shown. And we consider 
nonfactorizable spectator contribution in SCET, and find that they are also 
factorized to all orders in a s . The basic idea of the factorization properties 
in B decays into two light mesons was sketched in Ref. [16]. In this paper, 
we extend the argument and discuss intricate characteristics of SCET in non- 
leptonic decays, the details of the procedure of matching, and present all the 
Wilson coefficients and technical details in the calculation. We also present 
the analysis of B — > nn decays as an application, which is consistent with the 
calculation in the heavy quark mass limit at lowest order in a s . 

The organization of the paper is as follows: In Section 2, we consider the 
effective theories SCETi and SCETn, which we employ in the analysis of non- 
leptonic B decays and explain the field contents and discuss how to match 
these effective theories. In Section 3, we construct four-quark operators in 
SCET l5 which are gauge invariant by integrating out off-shell modes. This is 
achieved by attaching gluons to fermion legs and by integrating out off-shell 
intermediate states. We discuss the factorization properties of the operators. 
In Section 4, we match the full theory and SCETi and obtain the Wilson 
coefficients of the four-quark operators. In Section 5, we consider nonfactor- 
izable spectator interaction, in which subleading operators are enhanced to 
give the leading-order result. The structure of the subleading operators and 
the mechanism for the enhancement are determined by the power counting 
and the matching between effective theories. In Section 6, we consider the 
contributions to the heavy-to-light form factor. The time-ordered products of 
subleading operators from the heavy-to-light current and the interaction of 
an ultrasoft (usoft) quark and a collinear quark contribute to the form factor. 
In Section 7 we combine all the results to apply SCET to nonleptonic decays 
B — > 7T7T. And finally a conclusion is presented. In Appendix A, we show how 
the auxiliary field method can be used to derive gauge-invariant four-quark 
operators in SCETn- In Appendix B, we employ the auxiliary field method to 
derive gauge-invariant subleading operators in SCET!. 



2 Construction of the effective theories SCETj and SCET n 



The formalism of SCET and the procedure of the two-step matching are dis- 
cussed comprehensively in Refs. [9,15] and in Ref. [17]. We will not discuss 
them in detail, but we review them briefly here. The momentum of an energetic 
quark moving in the light-cone direction n M can be decomposed as 

m fit 1 

P» = ^n-p + pl + -n-p = O(\ ) + 0(A) + 0(A 2 ), (1) 
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where n**, n M are two light-like vectors satisfying n 2 = n 2 = 0, n ■ n = 2. 
There are three scales n ■ p ~ Q, pj_ ~ QX, and n ■ p ~ QA 2 , where Q is the 
hard scale. The momentum squared p 2 is typically of order QA, where A is a 
typical hadronic scale A ~ 0.5 GeV. We introduce a small expansion parameter 



A ~ yA/Q <C 1 to facilitate the power counting. We construct SCET in two 



steps, namely SCETi for y/QA < fj, < Q, and SCET n for // < \/C$A~. Formally 
we can construct SCET n directly from the full theory, but it is conceptually 
easy to construct SCETi and SCETn successively. 

In SCET, the collinear quarks £ and gluons A% have momenta p M = (n ■ 
p,n ■ p,p±) ~ Q(X 2 , 1, A). There are the soft fields q s , and A%, with momenta 
p^ ~ QX, and the usoft fields q us , A£ s with momenta p£ s ~ QX 2 . The collinear 
quarks interact with collinear gluons or usoft gluons, but not with soft gluons 
since they make collinear particles or heavy quarks off the mass shell. In order 
to derive the power counting method in SCETi, we move all the dependence 
on A to the interaction vertices, and determine the A dependence of all the 
fields [9]. Then we can construct all the operators in powers of A systematically 
by integrating out the off-shell modes of order Q. The guiding principles of 
constructing operators in SCET are gauge invariance [15], reparameterization 
invariance [13,14] and the power counting [18]. The Wilson coefficients of the 
operators can be obtained through the matching between the full theory and 
SCETj by requiring that the matrix elements of operators in both theories 
be the same at any order of the perturbation theory. There is no mixing of 
operators with different powers of A through radiative corrections as long as 
Q > A since the matching is performed perturbatively. 



In SCETn for fi < y/QA, we integrate out all the off-shell modes of order 
y/QA. Then the collinear fields in SCETn have momenta p^ ~ Q(X' 2 , 1, A'), 
and the soft degrees of freedom have momenta p s ~ QX', where A' ~ A/Q is 
a new small expansion parameter in SCETn- Here the collinear modes have 
momenta p 2 ~ A 2 , which is appropriate to describe hadrons of size ~ 1/A. 
The degrees of freedom of order ~ \/QA in SCETj are all integrated out, and 
the usoft degrees of freedom in SCETj of order p us ~ QX 2 remain in SCETn 
with momentum p s ~ QX 2 = QX', and we call them soft degrees of freedom 
in SCET n . The matching of SCE^ onto SCET n is performed at n ~ y/QA. 

A subtle technical point in matching is to require gauge invariance. In order 
to see how to keep the gauge invariance, let us consider a simple example of 
a heavy-to-light current. In the full theory the current qYb is matched to the 
operator in SCETi as 




C(n ■ V)£TWh, 



(2) 
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where C(n-V) is the Wilson coefficient, and = \n^n ■ V + V± is the label 
momentum operator. The factor W is the Wilson line which is given by 



w = ^2 ex P 



perm 



n ■ V 



(3) 



This Wilson line is obtained by attaching collinear gluons to the heavy quark 
and integrate out the off-shell modes of the intermediate states of order Q. 

In SCETn, the emission of soft gluons p^ ~ QX' makes the collinear fields 
off the mass shell and we have to integrate out the off-shell modes. This is 
achieved by attaching usoft gluons to external fermion lines in SCETj and 
by integrating out the off-shell modes of order QX 2 = QX'. It corresponds to 
factorizing the usoft-collinear interactions with the field redefinition [15] 

X 

e (0) = Y^, = Y^A n Y, Y(x) = Pexp(i0 J dsn ■ A us (ns)). (4) 

— oo 



The collinear effective Lagrangian can be written in terms of and A^ 11 , 
that is, the collinear-usoft interactions are factorized with the field redefini- 
tions given in Eq. (4). Then we match SCETt onto SCET n at a scale /i ~ y/QA. 
For example, the heavy-to-light current is matched as 

qTb -> C{n ■ V)lWTh -> C(n ■ V^W^TY^h, (5) 



where W (0) = Y^WY . In SCETn, the fields with the superscript (0) are fun- 
damental objects. We rename the usoft field as the soft field in SCETn. We 
write Y^h — > S^h, and lower the off-shellness of the collinear fields. 

Since the leading collinear Lagrangian is the same in SCETi and SCETn, we 

can simply replace C(n-V)^W^ — > C(n- V)^ U W U , and we obtain the final 
form of the operator 

qTb -> C(n ■ V)f I W Il TS^h, (6) 



where the superscript II, which denotes SCETn, will be dropped from now 
on. From the two-step matching, it is clear why the Wilson coefficients do not 
depend on the soft momentum. 

The advantage of the two-step matching is manifest when we consider time- 
ordered products in SCETj since these can induce jet functions involving the 
fluctuations of order p 2 ~ QA. In SCETj we can compute the jet functions 
with a well-defined set of Feynman rules independent of the computation of 
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the Wilson coefficients at p 2 = Q 2 . And the scaling of operators in SCETn 
is constrained by the power counting in SCETi and especially, SCETj puts 
a constraint of the number of factors of 1/A, which can be induced from the 
fluctuations l/QA. Therefore we can know how many powers of 1/QA appear 
for a given process at a given order in a s , and this power counting is not spoiled 
by the loop effects since there is no dependence on the soft momentum in the 
Wilson coefficients. This feature will be seen explicitly when we evaluate the 
decay amplitudes of nonleptonic decays which involve time-ordered products. 

In nonleptonic decays into two light mesons, there are collinear quarks moving 
in opposite directions. We choose these two directions as n M and n M . For an 
energetic quark moving in the n M direction, the momentum is decomposed as 

p£ = +PHL + = 0(A°) + 0(\) + 0{\ 2 ). (7) 

We denote a collinear quark by £ (x) which moves in the n M (n^) direction. 
These fields satisfy the relations 

« = 0, «{ = {, fe = 0, * X = * (8) 

The collinear gauge field in the (n M ) direction is written as (A£). The 
effective Lagrangian for x an d ^tc-, which can be obtained from the colllinear 
Lagrangian for £ and by replacing £ <-> Xi n ^ ^ respectively. This 
situation is simlar to two jets in the opposite direction [12], but we have a 
heavy quark interacting with collinear gluons in both directions, which makes 
the analysis more interesting and complicated. 



3 Construction of gauge-invariant four-quark operators in SCET 

The effective Hamiltonian for B decays in the full theory is given as 

tfe ff = % £ + c 2 oi + £ CiOt) , (9) 

V ^ p=u,c i=3,-,6,8 

where the local AB = 1 operators are given by 

Oi = (p a b a ) V ^ A (dpPi3)v-A, Of = (Pf3b a )v-A(d a pp) V - A , 

3 = (d a b a )v- A J2(qp<lp)v-A, A = (dpb a ) v ^ A Y^{q a qi3)v-A, 
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{d a b a ) V - A ^2{qpqp)v+A, Oq = {df3b a ) V - A ^2{q a qf3)v+A, 




'8 



■m b d a <T^(l + j 5 )(T a ) aP bpGl 



(10) 



Here p = u,c and d denotes down-type quarks, is the chromomagnetic 
field strength tensor, and T a are the color SU (3) generators. 

The process of obtaining the gauge-invariant operators in SCET requires two- 
step matching [19,20] because the SCET involves two different scales /i ~ m b 
and p,Q ~ \?m b k. First we match the full theory onto SCET: at /i — m b , and 
we match successively onto SCET n at /i . A concrete example of the two- 
step matching was illustrated in Ref. [21]. In order to construct the operators 
in SCET, we first have to specify which quark or antiquark goes to a certain 
direction to form a light meson. We set n M as the direction of a quark- antiquark 
pair to form a light meson, and as the direction of the remaining quark 
which combines with a spectator quark in a B meson to form another meson. 
Therefore the construction of the operators is process-dependent in the sense 
that we first specify the direction of each outgoing quark, and the number of 
operators in SCET is doubled because we have two possibilities to assign two 
quark fields in both directions. 

A generic four-quark operator for nonleptonic B decays in SCET has the form 
(£ri/i)(%r2x), or (;\Ti/i)(£r2x), where Ti and T 2 are Dirac matrices, and h is 
the heavy quark field in the heavy quark effective theory. These operators are 
derived from the operator q 1 Tib-q 2 T 2 q3 in the full theory where qi (i = 1, 2, 3) 
are light quarks. The operator {^Tih)(xT2X) * s obtained by replacing q\ by £, 
and q2 and q% by X- F° r the operator (x^ih)(^2X)j we replace q 2 by £ and q\ 
and q$ by the x fields. The second operator produces a light meson, in which 
one quark comes from the heavy-to-light current and another antiquark from 
the light-to-light current to form a meson in the direction. A remaining 
quark goes to the n M direction. We can transform this operator to the form 
(£r' 1 /i)(xr / 2 x) by Fierz transformation. In order to simplify the organization of 
the computation, all the operators will be written in the form (^Tih)(xT2X)- 
If a Fierz transformation is necessary to obtain this form, we apply Fierz 
transformation in the full theory and perform matching accordingly. Therefore 
we only have to consider operators of the form (£r 1 /i)(xr 2 x) ) but the two types 
of the operators with or without Fierz transformation are regarded as distinct 
and their Wilson coefficients are different. 

Though we know the form of the operators, the operator (^ih)(x^2x) itself is 
not gauge invariant under the collinear, soft and usoft gauge transformations. 
In order to obtain gauge-invariant operators in SCETn, we employ two-step 
matching. We consider the emissions of collinear gluons A£ and from ex- 
ternal fermions in the full theory and integrate out the off-shell intermedi- 
ate states of order ~ m b to obtain the collinear gauge-invariant operators in 
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Fig. 1. QCD diagrams attaching two gluons to external fermions to integrate out 
off-shell modes. The external gluons A n , An or A s make the intermediate states off 
the mass shell. Diagrams with gluons attached to other fermions are omitted. 

SCETj. Then we go down to SCET n , in which we consider the emissions of 
soft gluons which have momenta of order ~ m b \' = m b X 2 . The emission 
of soft gluons can cause intermediate states off the mass shell by the amount 
p 2 ~ m^A, which should be integrated out to obtain the final form of the 
operators in SCETn- To simplify the calculation, we will directly match the 
full-QCD operators to the operators in SCETn by considering the emissions 
of the collinear and the soft gluons from each fermion in the full QCD and 
integrate out all the intermediate states of order ~ nib and \fm h k to obtain 
the gauge-invariant operators in SCETn, as was done in Ref. [15]. The result is 
equivalent to the two-step matching described above. We will use the two-step 
matching explicitly in treating the time-ordered products. 



We can construct gauge-invariant operators by attaching collinear or soft glu- 
ons to each fermion in the full QCD and integrate out off-shell modes. Note 
that the soft gluons here have momenta of order ~ m^A' = m^A 2 , defined in 
SCETn- Because the ordering of gauge fields is important due to the non- 
abelian nature of the gauge fields, we consider corrections to order g 2 . Typical 
Feynman diagrams at order g 2 are shown in Fig. 1, and the remaining dia- 
grams in which gluons are attached to other fermions are omitted. But it is 
straightforward to attach two gluons to all the fermion lines making interme- 
diate states off-shell. 



If collinear and soft gluons are attached to the heavy quark [Fig. 1 (a), (b)], 
the intermediate heavy quark and gluon states are off the mass shell and we 
integrate them out. For the collinear quark £, the interaction with A^pi A^ 
makes the collinear quark off the mass shell. And the interaction with A^ or 
A^ makes the x field off the mass shell. We add all the possible configurations 
in which the intermediate states become off the mass shell and collect the 
terms at leading order in A. 

Let us introduce the factors 

A = B = C = D = (11) 

n ■ q n n- qn n ■ q s n ■ q s 
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where (q£) is the momentum of the collinear gluon (A£), and q% is the 
soft momentum of the soft gluon A%. The Wilson lines corresponding to each 
type of gluons are obtained by exponentiating the above factors as 



W= ^2 ex P 



perm 



s = J2 ex p 



perm 



n-V 
1 

n-TZ' 



n ■ A r , 
n- A, 



, w = J2 ex P 

perm 

> s = J2 ex p 



-9 



-9 



1 

-p^n ■ An 

n ■ Q 

i 

-n ■ A, 



perm 



n-K 



(12) 



Here = n ■ Vn^ /2 + (Q M = n ■ + is the label momentum 

operator for collinear fields in the n M (n M ) direction, and the operator 1Z is the 
operator extracting the soft momentum from soft fields. 

When we sum over all these diagrams, the color singlet operators of the form 
(L r i /l «)(x / 3 r 2X/3) and the nonsinglet operators of the form (^ /3 T 1 h a )(x a r 2X(s) 
are affected by the gauge fields differently and the final form of the four-quark 
operators can be written as 



8 = {yT x h a )%T 2X p) - H^L^F^ix^xs), 
On = (£pT 1 h a )(x a T2Xfi) -»• H^L^ s (^ a T l h p ){x 1 T 2 X5), 

where H°p, L° s (O = S, N) are the color factors. 



(13) 



One interesting case arises when we attach A^ and A^ to a heavy quark. At 
leading order in A, the amplitude in Fig. 1 (a), with A 1 ^ and A^ carrying the 
incoming momentum q£ and q£ respectively, is given by 



M a = 



m b (n- q n + n- qn) +n- q n n- qn 

IY _ iM\ n • A n 

xgri I m b n ■ An + n- q n n ■ An— — 

LV 4 ' n ■ q n 

(a -a W i \ n ' An} 

+ m b n ■ A n + n- q^n ■ A n —\ b, 

V A ' n. • fkr J 



4 ' n ■ qn ■ 

and the amplitude for Fig. 1 (b) with A% and A^ is written as 



(14) 



M b = % U^Tj^^b '' f, t+- 

2 n-qn-n-q n m b (n ■ q n + n ■ q^) + n ■ q n n ■ q^ 

xqT^n ■ A h ^^(m b + n ■ qjf)b. (15) 
n ■ q n v 4 > 

At first sight, these amplitudes contain complicated denominators which can- 
not be expressed in terms of A, B, C or D. However, if we add M a and M h , 
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we obtain 

M a + M b 



-7rfabcqTiT a = 6 + 5- gFi — 6, 

2 n • g w n • q n n ■ q n n ■ q^ 



(16) 



where there appear only A and B. Therefore the role of the triple gluon vertex 
is critical not only in determining the order of the Wilson lines but also in 
making the final expression simple. We will derive gauge-invariant operators 
using the auxiliary fields in Appendix A, and this cancellation is important in 
constructing the Lagrangian for the auxiliary fields. 

For singlet operators, when we sum over all the Feynman diagrams, we have 



H a(3 = 



g (-A + D)-g 2 AD L 



S 
7# 



5. 



76 • 



(17) 



to order g 2 . Here we show only the products of two different gauge fields, since 
they indicate the ordering of the Wilson lines. From Eq. (17), H^p suggests 
the Wilson lines to be (WS^) a p. For color nonsiglet operators, we have 



TjN 



7/3' 



g (-B + C)-g 2 BC 
g(-A + B-C + D) 

+g 2 (AC + DB -CB - DC - AB - AD) 



aS ' 



(18) 



where H^p suggests the Wilson line in the order (WS ) 7 p, and L^ s suggests 

the form (WS^SW ) a g. Therefore the gauge-invariant singlet and nonsinglet 
operators in SCET n are given by 



Os={(ZW)aT 1 (&h) a \ ■ [(xW)pT 2 (w\)p 
Os^WS^pY^h)^ ■ [{xW) a Y 2 {W\)p 



(19) 



The operator N can be written as ((£W S^pT^ h) a )((xW) a T 2 (SW 1 X )p) , 
but this is identical to N since (WS^S) ay <g> (PF f ) 7/3 = (WS ] ) ai <g> (SW^) y p. 

All the four-quark operators for B decays are of the form 0$ or On with dif- 
ferent Dirac structure. And the form of the operators O s and N in Eq. (19) 
manifestly shows the factorization of four-quark operators at leading order in 
SCET and to all orders in a s . In the operators Os and On, the interactions 
of A^ and A^ occur only in the heavy-to-light current sector, while the inter- 
actions of A^ occur only in the light-to-light current sector. If a collinear or 
a soft gluon is emitted from one sector and is absorbed by the other sector, 
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the interaction vanishes, or if it does not vanish, the momentum transfer is 
of order ^Jrr^K ) which is already integrated out to produce the Wilson coeffi- 
cients or jet functions. That is, the gluon exchange between the two current 
sectors is not allowed, and the gluon exchange is allowed only in each sector. 
Due to this property, the form of the operators is preserved even though there 
are any possible exchange of gluons to all orders. The Wilson coefficients may 
be different at different orders of the perturbation theory. This is an explicit 
proof of color transparency in SCET, and we can safely calculate the matrix 
elements of the operators in terms of a product of the matrix elements of the 
two currents. 

Note that the terminology "factorization" is used in two ways. First, it means 
that the matrix elements of the four-quark operators reduce to products of 
the matrix elements of the currents, and the matrix element can be written as 

(M.M^j^j^B) = {M 1 \ 3l \B){M 2 \j 2 \Q). (20) 



This was first assumed in naive factorization. It is possible only when there is 
no gluon exchange which connects the two currents and it is explicitly shown 
in SCET to all orders in a s . 

Another meaning of factorization appears in high-energy processes in which 
a physical amplitude can be separated into a short- distance part and a long- 
distance part. For example, exclusive hadronic form factors at momentum 
transfer Q 2 3> A 2 factor into nonperturbative light-cone wave functions for 
mesons, convoluted with a hard scattering kernel T as [22] 

/ fb f 

F(Q 2 ) = J dxdyT(x, y, fi)(f> a (x, fi)(/> h (y, fi) H . (21) 



We discuss both types of factorization in this paper. So far, we have considered 
the factorization of matrix elements. The second meaning of factorization will 
be discussed when we consider spectator contributions. 

By matching the four-quark operators in the full theory onto SCET, the ef- 
fective Hamiltonian for B decays in SCET can be written as 



G 

H c S = —7= E [K&K*d(ClT0iT + C 2 t0 2 t) 
V * T=R,C 

+ E v pb v; d £ cf T o iT 

p=u,c i=3,-,6 



(22) 



where 



11 



1R : 
2 R 

3R 
O m 
5R 



{rW) a {S^h) a ] v _ A \{^W)^w\ u 



V—A 1 



V _ A [(^W) a (W^) p 

(l d wusih) a ] v ^[(rwuw\<>) f: 



(fwsts)^)*] v ^ A £ [(x q wuw\% 
{l d wu^h) a ] v _ A Y\{rw) p {w\« 



V—A 1 
V-A> 



V—A 1 



V+A 1 



^ d wtfs) p (sh) a )] v _ A Y\(Tw) a {w\<>) 
\l d ws^s) p (s ] h) a } v ^ A [{x u w) a {W ] x u )p 



V+A 1 



\v-a 



o lc 

2 c=[(tw) a (S^h) a ] 

o 3C =Y;[(Vwsisus i h) 



V-A' 



V—A 1 



V-A 



(x d w%(^V) 



V-A 1 



V-A 



V-A 



(23) 



Here the summation over q goes over to light massless quarks. Since we specify 
the direction of each quark, we have the original operators iR , obtained from 
Eq. (9), in which the light-to-light current forms a meson in the n M direction. 
When a light quark in the heavy-to-light current moves in the n M direction 
and forms a meson with an antiquark in the light-to-light current, we use the 
Fierz transformation first in the full theory and match the operators. 

Note that there are no operators 5C and 6C - in SCET, which are obtained 
by Fierzing 5R and 6R . Neglecting color structure and the Wilson lines, 5C 
and Oqc have the form — 2£(1 + 75) h ■ x(l — 75)%, which is identically zero at 
leading order in SCET because 



X(l + 7s)x = + 7 5 )X = X(l + 75) ^X 



0, 



(24) 



since ifix = 0. In the literature [7], the effects of the operators O^c and 0%c are 
sometimes known as chirally-enhanced contributions. Even though the effect 
is formally suppressed, the numerical values may not be negligible. When the 
equation of motion for the currents is applied, there is an enhancement factor 
m< M / m g m b, where m M is a meson mass and m q is the current quark mass. But 
in SCET there is simply no such operator as 5 c and Oqc at leading order. 
This contribution is formally suppressed in powers of A in SCET, but the 
coefficient of these operators can be large. In phenomenological applications, 
we have to know how to treat the chirally-enhanced contributions, but we will 
not consider them here. 
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4 Matching and the Wilson coefficients in SCET 



The Wilson coefficients of the four-quark operators in SCET can be deter- 
mined by matching the full theory onto SCET. We require that the matrix 
elements of an operator in the full theory be equal to the matrix elements 
of the corresponding operator in SCET. We use the MS scheme with naive 
dimensional regularization scheme and anticommuting 75. All the external 
particles are on the mass shell. 

The radiative corrections for the four-quark operators in the full theory are 
shown in Fig. 2. As a specific example, the amplitudes for Fig. 2 (a) to (d) for 
the operator 0\ = (du)v-A(ub)v-A are given as 



= ^ Ta ^k( T a)il(djUk)v-A(Uibi) V ^ A 



X 



^uv 



12 xim b 
- + -(In 

e z e fx 



2 — 3Xi 7T 

f lnxi - 2Li 2 (l - x{) - — 



A _ 4 + 21n^-21n 

> fx 

2 



2 x x m h 
fx 



1 — Xi 



12 

iM b (1) = ^{T^ jk {T a ) a i^jU k )v-A{uibi)v-A 



x 



2 



+ i-'(ln^-l)-5 + 41n^ + 21n 2 ^ 
e 2 e v jx ' jx fx 



7Y 



— lnx 2 + 2Li 2 (l-x 2 ) + — 
x 2 12 



(a) b 




(b) 




(c) 




(d) 




(e) 




(f) 




(g) 



(h) — #o7Jolrl — 



Fig. 2. Radiative corrections at one loop in the full theory. The momenta pi, p 2 , p 
are outgoing with Pb = P + Pi + P2- Infrared divergences exist in diagrams (a)-(f). 
Diagrams (g) and (h) are infrared finite. In (h), the square is the chromomagnetic 
operator Og. 
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iM c (1) = ^(T a ) jk (T a ) u (djU k )v-A(v4bi)v-A 



x 



i 2 2/ —xximl\ _ , nz-np^ 

+ 3 + 7( 2 - ln ^r J; )- 1 + ln (^r^) 



euv 



/i 



7T 

~6 



= ^-(r a ) jfe (r a )j/(d j u fe )y_ j4 (^&i)v-A 

47T 



X 



1 2 2/, 
---(2 -In 

euv e e v 



-xx 2 m£\ 



-8 + 31n -^ -ln 2 -^ + ^ 



6 



(25) 



where Li 2 (x) is the dilogarithmic function. And x, Xi j2 are the energy fractions 
given by x = n-p/rrib, Xi = n ■ Pi/rrib for i — 1, 2. The prescription for ln(— Xj) 
is given by ln(— Xj — ie) = lnxj — i7r. If we add all these "nonfactorizable" 
contributions, the infrared divergence cancels and the only infrared divergence 
comes from the vertex corrections of the currents [Fig. 2 (e) and (f)]. Since 
the vertex correction of the light-to-light current is the same both in the full 
theory and in SCET, it cancels in matching. And Fig. 2 (f) is given by 



lM f ] = ^ L (du) v _ A (ub) v _ A 



X 



+ 



An 
1 



x 



X 



1 2 / xrrib 
- + -(In 



lnx — 2Li 2 (l — x) 







3 + 21n^ 



2 In 



2 xrrib 



71 

12 



(26) 



For O5 = (^jy-iEgtely+A) the corresponding amplitudes are given as 



iMf = ^(r a ) J - fc (r a ) a (d i 6 I )v-A(g j g fc )v + A 



X 



iM, 



(5) 



+ 



ct 



euv « 
2 



1 2 /, 
^ + 7( ln 



XiTTlb 



Xi 



e x /i 
lnxi — 2Li 2 (l — Xi 



l)+2-41n^-21n 

2 



2 xim b 
fx 



7T 
12 



^ (T a )jk {T a ) ii ( dibi ) y_ A (<?j <?fc ) U+A 



X 



+ - 

euv e 
2-3x 2 



x 2 m b 
/i 



2 x 2 m b 



l)+l_21n^ + 21n 2 
' /i /i 

J2_ 



71 



lnx 2 + 2Li 2 (l - x 2 ) + 
iMf = -^(T^kiTaUdMh-Aiqjqkh+A 
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X 



1 2 2/, 
_+ + (2-ln- 



-xx\mi 



+ 5-31n-^ 



(5) 



+ ln 



2 -ra^j 



7T 2 1 



47T 



/i 2 6 
(T a )j k (T a )u(dibi) V -A(qjqk)v+A 



x 



4 2 2 /, 
( 2 -ln 

euv e e v 



-xx 2 ml\ 



-2-ln 2 



-XX 2 ml 7T 2 

+ ~6~ 



(27) 



and iMj 5 - 1 is the same as iM^p except the Dirac structure. The infrared diver- 
gence of the nonfactorizabie contributions in Eq. (27) cancels, and the infrared 
divergence from the vertex corrections is cancelled in matching. 



Using the color identity 
1 



(T a )jk(T a )u 



8jl5ik — 2j-j:8jk$il, 



(28) 



the operators in Eq. (25) becomes 0\j1 — 2 /(2N), and the operators in 
Eq. (27) becomes 5 /2 — 6 /(2N). In the radiative corrections of 2 , the 
operator becomes C F 2 due to the color factors in Fig. 2(a) and (d), while 
it becomes Oi/2 — 2 /(2N) from Fig. 2(b) and (c). Therefore the radiative 
corrections for 2 , or for nonsinglet operators in general, have infrared diver- 
gences in all the diagrams. 

The corresponding radiative corrections for the four-quark operators in SCET 
are shown in Fig. 3. For singlet operators, the radiative corrections exist only 
in the heavy-to-light current sector [Fig. 3 (a), (b)], in which the infrared 



x 





Fig. 3. Radiative corrections at one loop in SCET. Curly lines with a line represent 
collincar gluons, and curly lines represent soft gluons. 
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divergence appearing in the calculation is exactly the same as the infrared 
divergence in the full theory. The diagram (c) is included, but since the result 
is the same in the full theory, it is cancelled in matching. For nonsinglet 
operators, we need all the diagrams in Fig. 3 with additional collinear gluon 
interactions in the light-to-light current sector [Fig. 3 (e) and (f)] and soft 
gluon interactions [Fig. 3 (b) and (d)]. When we calculate these diagrams, the 
infrared divergence is exactly the same as that in the full theory. Therefore 
we can safely match the full theory onto SCET since the infrared divergence 
cancels, and the Wilson coefficients can be calculated. 

In matching the theory at one loop, we calculate perturbative matrix elements 
in the full and effective theories. All the long-distance physics is reproduced 
in the effective theory, and the difference between the two calculations deter- 
mines the short-distance Wilson coefficients. We treat both ultraviolet and 
infrared divergences using dimensional regularization, with the final collinear 
quark on the mass shell. In this case all the Feynman diagrams in SCET are 
proportional to 1/euv — l/ e m — 0. The ultraviolet divergences are cancelled 
by the counterterms in SCET, and all the infrared divergences cancel in the 
difference between the full theory and SCET. Therefore the Wilson coeffi- 
cients of various operators in SCET can be obtained by calculating radiative 
corrections in the full theory. In the full theory we also have to consider the 
Feynman diagrams with fermion loops and the effect of the chromomagnetic 
operator. These contributions are included in Fig. 2 (g) and (h). 

The Wilson coefficients are, in general, functions of the operators n ■ V/m b , 
n- Q/rrib and n ■ Q) /m b . The matrix elements of the four-quark operators can 
be written in terms of a convolution as 



/— / n ■ \ 
dxdx x dx 2 C{x, x u x 2 ){£W6[x - —^—jF^h 



xY W5( Xl - !^)r 2 5(x 2 + — WV> (29) 

Due to the momentum conservation x, x\ and x 2 satisfiy x + X\ + x 2 = 2. For 
nonleptonic decays into two light mesons at leading order in SCET, we can set 
x — 1, x± — u, and x 2 —u = 1 — Xi, and the matrix element can be written as 

M —> J dr] CWgWr^h-xWSiri - Q + )T 2 W\). (30) 



where u = X\ = 1 — x 2 = n/ (4i?) + 1/2 with Q + — n ■ Q} + n ■ Q. However, we 
list all the Wilson coefficients for general x, x\ and x 2 , which will be useful for 
other decay modes, or nonleptonic B decays at subleading order in SCET. 
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By adding the wave function renormalization of the heavy quark, the Wilson 
coefficients at next-to-leading order are given as 



C\R — 


1 - 


M 

2N 


+ C F A 2 


Ci 


+ 


^3 

— L/ 2 , 2 ft 




+ 


1 - 


r<p — 


1 - 


A 1 
2N 


+ C F A 2 


c 3 


+ 


A* 

— 4 , L/ 4R 


2 3 


+ 


1 - 


riP — 
L/ 5R — 


1 - 


A 5 
2N 


+ C F A 2 


c 5 


+ 




2 5 


+ 


1 - 


Cic = 


1 - 


A 3 
2N 


+ C F A 4 


Ci 


+ 


-^-C 2 , C 2 c 


_A 3 

~Y Cl 


+ 


1 - 


nP — 
°3C — 


1 - 


A 3 
2N 


+ C F A 4 


C 3 


+ 


Ai r 1 
2 4_ 2iV 


Cf, 







M 

2N 
As_ 

2N 
A_ 

2N 

A 

2N 



A* 

r<p — 3 r 1 _i_ 
°4c — ^~°3 + 



^-^ + C F A 2 ]C 4 + 1 -Cl 



+ C F A 4 C4, 
+ C F Aj Cq, 

c 2 

(31) 



where Cj's are the Wilson coefficients from the full theory. And the coefficients 
Ai to order a s , evaluated at ji = m b are given as 



A 2 : 

Ay- 



A, 



A 1 = -± -18 + 3m(-x) +ln 

4:71 L 



X 



X\ 2 — 3xi 



In— + 



X\X 2 X 2 1 — X\ 



lnxi + 



1 - 3x 2 
1 - x 2 



In x 2 



-2Li 2 (l-xi)+2Li 2 (l-x 2 ) 



47T 
47T 



2 — X 7T 2 

-5-21n 2 x + In a; - 2Li 2 (l - x) - — 

-9 + ^— - lux - 2Li 2 (l - x) + 2 In 2 x 2 + 2 ln 2 (-xi) - In 2 — — 
1 — x x 



2 7T 2 

In x 2 + 2Li 2 (l - x 2 ) - — 

l-x 2 6 



An 



2 — 2>X\ 

— 14 — 2 In 2 xi — ln 2 (— xx 2 ) + 3 ln(— xx 2 ) + lnxi 

J- v6 1 



7T 



-2Li 2 (l-x 1 ) + - 



47T 



6 — 3 ln(— x) + In 



x 



In 



1 - 3xi . 2 - 3x 2 

in xi m x 2 



2Li 2 (l-x 1 )+2Li 2 (l-x 2 ) 
3 + 21n 2 x 2 + 21n 2 (-Xi) - In 2 



An 
2 - 3x 2 



XiX 2 X 2 1 - Xi ' 1 - x 2 

xi 1 — 2x 



X 



l-x 



lnx — 3 ln(— Xi) 



lnx 2 -2Li 2 (l-x)+2Li 2 (l-x 2 )-- . 
1 - x 2 6 J 

A 7 = ^[-2 - 21n 2 xi - ln 2 (-xx 2 ) + - — — lnxi 
An 1 1 — Xi 
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-2Li 2 (l-a:i) + 



7T 

12 



(32) 



And the contribution Cf from fermion loops and the chromomagnetic operator 
[Fig. 2 (g) and (h)] at /i = m b are given by 



nP _ &s r 



'C7i(^-G(a p ))+C7s(^-G(0)-G(l 



+C 4 (^-3C(0)-C(s c )-C(l)) 
+C 6 (-3G(0) - G(s c ) - G(l)) - (C 5 + C 8 ) — 

where s p = ml/ ml (s c = m\/ ml), and G(s) is given by 



Xi 



(33) 



_L 

G(s) = -A J dzz{\ - z)(s - z(l - z)(l - x x ) - 



(34) 



We can add C?/(2N) in Cf and C 5 R , and Cf/2 in and Cf with x x -> x 
in Eq. (33). But in physical processes in which the final-state consists of color 
singlet mesons, the contribution of Cf cancels, hence omitted in Eq. (31). 



5 Nonfactorizable spectator contributions 



In nonleptonic B decays, we also have to consider the spectator quark which 
can interact with the b quark or other quarks forming light mesons. In this 
section, we concentrate on the spectator quark interacting with the quark 
and antiquark pair (x fields) produced by the weak current, which we call 
nonfactorizable spectator contributions. The spectator quark can interact with 
the heavy-to-light current, which will be treated in the next section. Here 
we apply the two-step matching explicitly. The time-ordered products are 
constructed in SCETi, and they are evolved down to SCETn- 

Nonfactorizable spectator contributions arise from the interactions of the glu- 
ons A% in the light-to-light current with a spectator quark, which becomes a 
collinear quark as a result. However the operators Os and at leading order 
in SCET do not involve the interaction of Afi. Therefore we take into account 
subleading operators which involve A 1 ^ in the light-to-light current. And the 
Lagrangian describing the interaction of collinear and usoft quarks begins with 
0{\) compared to the leading collinear Lagrangian. But the propagator of the 
exchanged gluon is of order A~ 2 . Therefore the nonfactorizable spectator con- 
tribution is of the same order as the leading contributions from the four-quark 
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Fig. 4. A Feynman diagram for nonfactorizable spectator contributions from the 
subleading operator in the light-to-light current. The soft momentum k is incoming, 
and pi (i = 1,2,3,4) are outgoing. 

operators. In order to evaluate decay amplitudes at leading order in SCET, 
we need to include the nonfactorizable spectator contributions, as shown in 
Fig. 4. Here we consider a collinear gluon from the light-to-light current 
interacting with a spectator quark to produce a collinear quark £. 

The Lagrangian for the collinear and usoft quarks at order A is given by [21,23] 
4? = \ W ^nW] + h.c. (35) 



In SCETi, there are two independent subleading operators suppressed by A, 
in which collinear gluons Al^ come from the current X^2X- They are given by 



<>f> = ((^^^^^^^([Wt^^^j^^^x)^, (36) 

where the index % runs over all the possible forms of the operators as shown 
in Eq. (23). The Wilson coefficients of these operators are 1 at lowest order 
in a s . The form of the operators of~ a ' lb ^ can be obtained in a straightforward 
manner, but the ordering of the Wilson lines is nontrivial. It will be explained 
using the auxiliary field method in Appendix B. 

In SCETi, the nonfactorizable spectator contribution is given by the matrix 
elements of the time-ordered product 



r 2l (vy t x)/ 3 



1 — t. 



TP = I d 4 xT 



O^\0) + or\0),i£$(x) 



(16)/ 



.(1)/ 



(37) 



In order to go down to SCET n , we decouple the collinear-usoft interaction 
using the field redefinitions [15] 
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£(°) = yt£ A^ = Y^A n Y, Y(x) = Pex V (ig J dsn ■ A us (ns)) , 

— oo 

X 

X m = Y ] X , 4 0) =F%F, Y(x) = Pexp(tg J dsn ■ A us {ns)) . (38) 

— oo 

The collinear effective Lagrangian can be written in terms of x^ an d 
A^^, A^" 1 , that is, the collinear-usoft interactions are factorized with the 
field redefinitions given in Eq. (38). 

Matching at /i ~ ^/m^A, the small expansion parameter changes from A ~ 
^A/m b to A' ~ A/mb with a definite power counting procedure for any oper- 
ators. Recall that we rename the usoft field of order p ~ m b \ 2 ~ m b \' as the 
soft field in SCETn. In that sense the usoft fields become soft and the Wilson 
line Y becomes the Wilson line with soft gluons in SCETn but "soft" means 
the momentum of order m b \' . Without introducing addtional Wilson lines, we 
denote this Wilson line as S and in the matching we replace Y — > S without 
any ambiguity. And the operators are matched onto the operators in SCET n . 

The nonfactorizable spectator contribution comes from the matrix elements 
of six-quark operators. In calculating the matrix elements of we first 
factorize the usoft-collinear interactions using Eq. (38) with Y — > S, Y — > S, 
and project out color indices in such a way that the quark bilinears W^) 
and (xW, W^x) forming mesons are color singlets. Since the final expression 
involving the Wilson lines are nontrivial, we show explicitly how the color 
projection is performed for 0\ la \ At order a s , we will extract a collinear gluon 
from the factor (WUp n ±W) in C^ q and to contract them. Neglecting 

the Dirac structure and keeping color-dependent parts only, the time-ordered 
product of 0- 1<1 ' ) with C^ q has the form 



[(eW^)/j • ha] [(xWS^ST a S^) a ■ [SW ] X )p] [{q s SUT a ) l6 {W^)s] 

= [(W^)^ ■ h a ] [{xW) v (s ] ST a s^) ua ■ {S) Pp {W\) p ] 

x[(q s SUT a ) l5 (W^) 5 _ 



x(s 1 ST a S^JS^S) pp (T a ), s 

[(&?)(, ■ K] [(xW) u ■ (W\) v ] [faS) 7 • (W*fl/j] 



1 



C_ 

N 2 



7/( 



-§ [(tW) fi ■ (&h) a ] [{xW) v ■ [W\) v ] [{q s S) a ■ (W^)p] , (39) 
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where the dots denote the Dirac structure and the color projection is per- 
formed after the arrow in Eq. (39). The ime-ordered product with can 
be projected in the same way. The time-ordered products of the operators 
Of a ^ with the color singlet structure vanish due to the color structure. 

The matrix element of the time-ordered product (T- 1 ^) at order a s is given by 



' ! ' v = — A7ia s ^r I (111 ■ .V 



dn ■ k 
Air 



in-kn-x/2 



x{([(^)^r u (^) a ] [(xw),— 



n- Qt2 



x 



(l s S)a(n-x) 



1 1 ' iw^uo)]) 



n-TV '"n-V 



X 



1 



1 



n 



(40) 



Let us consider in detail the matrix elements in Eq. (40) for different Dirac 
structure Tu <g> r 2i . For simplicity we omit all the Wilson lines and the mo- 
mentum operators in the following calculation. For 7^(1 — 75) 07^(1 — 75), we 
can evaluate the first term in the curly bracket in Eq. (40) as 



(IM 1 - ■ x^ilYii - 7 5 )x • la-rib) 

= (^(1 - T*)ha ■ X^9T - tK)(1 - 7 5 )X • H^Zp) 
= 2(^(1 - 7 5 )h a ■ q a ^p ■ x\ (1 - lM 

-(^(1 " ls)K • x\llll (1 - 7s)x " q a 7ib)- (41) 

In the first line we can replace 7^ by 7^ and we use the Fierz transformation 
to arrive at the last relation. The first term in the last relation in Eq. (41) can 
be further simplified as 



2^(1 " ■ q al fa ■ x>-{\ - 75 )X 

= ^7^(1 -ls)K -q a [l^ -75) +7m( 1 + 75)]^/3 

= ^?7 M (1 - 7s)£ • 27^(1 - ls)h ■ xfi(l - 7s)x 



1 - 7s)x 
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= 1^(1 - 75 )e • ¥(1 - 7 5 )/i • - 7s)x- (42) 

The part proportional to 7 M (1 + 75) vanishes when we apply the Fierz trans- 
formation. In the third line we use the Fierz transformation for the product 
of the first two currents. Similarly, the second term in Eq. (41) is simplified as 



-f/j7«/(l - ls)h a ■ X27±7±(l - 7s)x • 5«7^/9 

4 

= -iM l - 75)^ • X^7"(l + 75)7^X • ^7^/3 

_ # 

= ^/37 M (1 - 7s)X 7 • 9a7Ai(l - 75)^/3 • X 7 ~(l - 7s)^a 



1 



W - 7s)£ • - 7s)x 7 • X 7 ( x + 7s) 



= --^${1 - 75 )C • qi{l - 75)/* • xfi(l ~ 7s)x, (43) 

where the Fierz transformations are applied successively. When we add Eqs. (42) 
and (43), they cancel. The second term in Eq. (40) is given by Eq. (43) with 
an opposite sign. As a result, for Yu <g> Y 2 i = j v (l — 75) <8> 7"(1 — 75), only the 
second term in Eq. (40) contributes and the matrix element is given as 



x(xW-^IFV(1-75)x>- (44) 

Note that the matrix element of the six-quark operators in Eq. (44) is written 
in a factorized form. It is because a collinear or a soft gluon in one sector 
cannot interact with the particles in the other sector. That type of interaction 
makes the intermediate states off the mass shell by ~ rrib or y/rribA. These 
off-shell modes are already taken care of in the Wilson coefficients or in the 
jet functions. Therefore the six-quark operators describing the spectator in- 
teractions are factorized and the matrix elements of the six-quark operators 
can be computed by the products of the matrix elements of each current. 

For 7j,(1 —75) 07^(1 + 75), only the first term in Eq. (40) contributes and the 
matrix element is given by 
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x(xW^^W^(l + l5 ) X ). (45) 
And for (1 — 75) (g) (1 + 75), the matrix element is given as 

hlwJ—W^{\ + l5 )0(q s S^—SW±(l - ls)h) 
16 n ■ P n ■ IV 

X(xW (^Qt - ^q)^V7^(1 + 7 5 )x>. (46) 

As can be clearly seen in the final expressions in Eqs. (44)-(46), the gluons A^, 
A£, and can be exchanged only inside each meson. This is true at higher 
orders of a s . Though the Wilson coefficients can be different, the structure 
of the operator is the same to all orders in a s . Therefore the nonfactorizable 
spectator contribution is factorized to all orders in a s . 

Let us calculate the matrix element t/ 1 ^ explicitly for Tu Cg>r 2 j = 7^(1 —75) <E> 
7^(1 — 75). It is given by 



C F r, rdn-k e ^-kn-x/2 



(Th = -4>™ aJ £ J dn ■ x J 



47r n ■ kn ■ p 2 n ■ p 3 



x l -{M^Wf{\ - 75 )^| )(M 2 |x^(l - 75)^X10) 
x(0\q s S(n-x)i(l- l5 )S^h(0)\B), (47) 

where we integrate out n ■ x and x± explicitly. The matrix element involving 
the B meson can be calculated as 



(0\q s S{n ■ x)i(l - ^)S ] h\B) = J dr + e~ tr + n - x/2 Tr [* B (r+)yi(l - 75) 

= Jf«p- J dr + e-^Tr [I^ 7 ^(i - 75 )]^( r+ ) 

= -if B m B J dr + e~ ir ^ x ' 2 <P + B {r + ), (48) 

where the leading-twist B meson light-cone wave function is defined through 
the projection of the B meson as [26,27] 

*b{t + ) = _^£[ili(^+(r + ) +^i(r + )) 75 ]. (49) 
And the light-cone wave function for the light mesons is defined as 
(M 2 \xWi l5 5( V - Q + )W\)\0) 
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1 

-ifu^E J du5[ V - (4u - 2)E]<f )M2 {u). (50) 



For 7^(1 — 75) ® 7 M (1 =F 75), the matrix element which is given by 



(Tj (1) ) — ^ 2 s /b/mi/m2™b 

x /dr + ^±> /^^M [dv^l. (51) 
7 r + J u J v 

For (1 — 75) (g) (1 + 75), the matrix element is zero if we use the leading-twist 
B meson wave function because 



1^(1 + -K)h\B) = _!^p tr + 75 )!±^0+ + ^-) 75 ; 



4 

if Bins 
8 



<^tr fflfrU = 0. (52) 



If we use higher-twist wave function for the B meson, there may be nonzero 
contributions, but this is expected to be suppressed. 



6 Spectator contribution to the form factor 



In Section 5, we have considered the nonfactorizable spectator contributions 
arising from the subleading operators in which we include only the sublead- 
ing part from the light-to-light current. However, there is also a subleading 
operator coming from the heavy-to-light current, but this contributes to the 
heavy-to-light form factor. It is considered first in Ref. [20], and we discuss in 
detail here in the context of nonleptonic B decays. 

The operators in SCETx, which contribute to the form factor, are given by 

4 0) = Gwr u h)(xWr 2i w*x), 

4 16) = (W(WU^ n± W)^-h) {xWT 2l W ] x). (53) 

Here we list only singlet operators. The nonsinglet operators give the same 
matrix elements as the singlet operators with the color suppression factor 1/N. 
We will consider only singlet operators from now on. 
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Fig. 5. Tree-level graphs in SCETi for the spectator contribution to the 
heavy-to- light form factor. The first diagram contributes to ^1^,4, and the second 
diagram contributes to ^1^4^^. 

The interaction of collinear and usoft quarks is given by [20] 

igl^--$ r lWq us + h.c, 4f = igl—^—$Wq us + h.c, 
in • D n in • D n 

i 4 i ^w^ mqus+h - c -> (54) 

= [in • Am 0> n ± ], 19$ = [in ■ D n , ip us + ^gn ■ A n \. (55) 

At leading order in SCET, the relevant time-ordered products are given as 

T F = J d 4 xT[4°\0)iCf q \x)}, T[ t = J d'xT[4 la \0)iC^(x)}, 
T£ = J d 4 xT[4 lb) (0)i££(x)}, T( t = J d'xT[4°\0)ijCf q b \x)}, 

T» F = J d'xT[4°\0)iLf q a \x)}, 

T» F = J d 4 xd 4 yT[4°\0)iC^(x)iCf q \y)}, 

T 6 N / = J d*xd*yT[jP(0)iC%(x)i££ q \y)}, (56) 

where £^ is the leading collinear Lagrangian and C^J is the subleading gluon 
Lagrangian, which can be found in Refs. [20,21]. The Feynman diagrams at 
lowest order in a s for the time-ordered products are shown in Fig. 5 schemat- 
ically. Compared to the case of the nonfactorizable spectator interactions, the 
leading-order heavy-to-light current jf^ contains A%, therefore it contributes 
to the heavy-to-light form factor starting from the leading order. 

As suggested in Ref. [20], we absorb the nonfactorizable parts T^ F (k = 4, 5, 6) 
into the definition of the soft nonperturbative effects for the form factors at 
this order. Among the factorizable contributions T t f , only T£ is nonzero at 



,(26) 



where 

igfl 
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order a s . For <g> T 2 i = 7^(1 — 75) <8> 7^(1 — 75), the matrix element of T 2 F { is 
given by 



(T£) = (M 1 \i J d 4 xT[jl lb \0) : jC,W(x)}\B) 

92 1 ■ (dn-x /*L* 
J J n ■ k 



47r rribn ■ p\ 

a s A-k 2 C f ., 2E r 0Mi(«) f dr A 



x ^( W(l - 75)^ • W(l - 75)W^X • g s S(n • x)i{\ - lb )S^h) 



iJmJm2Jb — m B / du / <P B {r+), ( 5 <) 



4tt JV mb J u J r + 

where we use Fierz transformations successively. For <g> Y 2 i = 7^(1 — 75) <8> 
7"(1 + 75), we obtain the same result as in Eq. (57). For (1 + 75) <g> (1 — 75), 
it vanishes. 

The form factors for B decays into light pseudoscalar mesons are defined as 



(P(p)\qYb\B(p B )) = f + (q 2 ) 



U „ m B ~ m P n 

Pb + P^- 2 V 



+/o(g 2 )^-^^, (58) 



where q 11 = p B — p^. In SCET, the form factor at order a s is given by 



C F hnJmiB 2E r q s (/i) 
N AE 2 m b J X r+ xr + 
+ (l + K)C(/i ,/i), (59) 



, fn , C F f M1 f B m B 2E r ajjj) + 

/ + (0) =ira s — — — — J dxdr + ——(j) M1 {x^)(j)l{r +) ^) 



where ((fiQ, fi) is a nonperturbative function introduced in Ref. [24] in large- 
energy effective theory [25]. A similar nonperturbative function is introduced 
in Refs. [9,13], based on SCET. The procedure in obtaining independent non- 
perturabtive functions in the form factor is different in the large-energy ef- 
fective theory and in SCET, but the number of independent nonperturbative 
functions is the same. And K at /x = m& is given by 

«=-^r(?+u)- (60 » 



The expression for / + coincides with the result in Ref. [20] with m B = 2E at 
order a S) and K is calculated in Refs. [9,13]. 
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7 Application to B — > txtx decays 



As an application, we consider the decay amplitudes for B — > mi, which can 
be written as 



(mi\H cS \B) = ^= £ V pb V,; d (nn\A p \B) , 

V ^ p=u,c 



(61) 



where the operators A p are given by 



v-4.p Oj 



{Ch) v ^x u )v 



(Ch) v -A(rx q )v- 



+ a\ 



(f^-A(xV)v-. 



(62) 



The notation a\[0] means that af are the sum of the amplitudes initiated by 
the operator O. They include the contribution from the operator itself, and 
the spectator contributions. That is, af [O] can be written as 



am=T[ + Nf + F[, 



(63) 



where Tf is the contribution from the four-quark operators, Nf is the nonfac- 
torizable spectator contribution, and Ff is the spectator contribution for the 
heavy-to-light form factor with a{ = a\ = 0. 

Before we present the decay amplitudes for B — > mi at order a s explicitly, we 
show the amplitudes af in general to all orders in a s , which are derived from 
SCET. The form of Tf from Eq. (30) can be written as 



Tf = J drj Cj Sji (rf,fi ,fx) 

x (£W ltM (l - 75)^/1 • xW5( V - Q + ) 7 M (1 T 1^0 x) 

= ± i f M 22E J ^Cj ffii K// o ,//)0M2K/i)(M 1 |^|(l-75)-5 t /i|S) 


1 

= ±im 2 B f M 2 J du((fx , fi)Cj Sii (u, fx , fj)4>M2{u, fx), (64) 


where the upper (lower) sign corresponds to % — 1, • • • , 4 {% — 5). The effective 
Wilson coefficients are evaluated at /i = m& and evolved down to fx = 
fio = y/rribA to be matched onto SCET n . Then they evolve down to the scale 
fx, where the matrix elements are evaluated. 
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The matrix element of the heavy-to-light current is given by [24,13,9] 

{M^Wfil - l5 )S^h\B) = 2m B (. (65) 

The nonperturbative parameter ( is matched at \i = /i and it evolves down 
to /i, and the wave functions are evaluated at ji. Note that we do not include 
the radiative corrections in (, since they are taken into account either in the 
Wilson coefficients or in the spectator contributions. 

The spectator contribution Nf involves the time-ordered products, which can 
be written as 

T[Ol la \x) + oi lb \x),iC^(0)]=5(^)5 2 (x ± ) J drjdr] J dr + e ir ^ x/2 

x J?(v,V,r+)O i (r},r},r + ), (66) 

where Jf are the jet functions, which are obtained by matching SCETi onto 
SCET n . The index i denotes all the possible operators for the decay. And the 
operator Oi is given by 

Oi(v, % r + ) = WruHv - v+)w^ ■ xw^siv - Q + )w ] x 

xq s i(l-K)$(rt-r + )tfh. (67) 
Therefore Nf can be written as 

N? = J d*xC^T[0t\x) + 0l lb \x),i£Mm 

= J dvdvdr+C^(no, v)J*(u, v, r+, fi , n)Nif B fMifM 2 

X0 M i(M, / u)0 M2 (v,/i)0B( r + ) / U ) ? (68) 

where i are the effective Wilson coefficients, and the fi dependence is shown 
explicitly Here iVj are the normalization constants when we evaluate the ma- 
trix elements of Oi. Since n-pMi = n-pM2 = 2E, the variables u and v satisfy 
the relations 

u = — + -, v = — + -. (69) 
AE 2' AE 2 1 ; 

The factorizable spectator contribution to the form factor can be written sim- 
ilarly as 

Ff = J dudvdr + C^ Sii (/j, , fi)J?(u, v, r+, /x , jj)Nif B f M ifM2 



28 



X(f>Ml(u, (J,)(j>M2(v, /i)0s(r + , /i), 



(70) 



where C^ Si are the effective Wilson coefficients, and Jf are the jet functions 
obtained from matching SCETi onto SCET n . The spectator interactions Nf 
in Eq. (68) and Ff in Eq. (70) are factorized into the short-distance and the 
long-distance parts to all orders in a s , and the convolution integrals are finite. 

Now let us calculate the decay amplitudes for B — > nit at order a s , based on 
the general expressions on Tf Nf and Ff. The contributions Tf are obtained 
by the convolutions of the following effective Wilson coefficients 




with the hadronic matrix elements of the four-quark operators. For the decay 
amplitudes at leading order, we can put x = 1 and let u = x\ and u = x 2 = 
1 — u, and the amplitudes Tf s, evaluated at /i — mi, are given as 



Tf^ m lCU[(C 1 + ^) { l + K) + ^p 2 F 
Tf^ m lCU[(C 2 + ^)(l + K) + ^ Cl F 



N ■ 

C; 



4tt N 



Tf = im 2 B CU\ [(C* + ^)(1 + K) + ~[CsF + C 1 (-- G(s p 



4tt N 



+ C 3 (^- G(0) - G(l)) + - 3G(0) - G(a c ) - G(l 

+ C 6 (-3G(0) - G(s c ) - G(l)) + G^ 8 (G 5 + C 8 )] |, 
T? = -imlCA {(C 5 + ^)(1 + K) + ^—C^-F - 12) , 
where C F = (N 2 - 1)/(2N), N = 3, n f = 5. And F is given as 



(72) 



F = -18 + fl = J dug(u)</> w (u), G wfi = J du^—^>M, 



q(u) = 3 In u — 3iir 

v 1 l-u 



1 — 3u 

2Li 2 (l — u) H In u + In 2 u — (u <-> u) 

1 ?x 



(73) 
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The nonfactorizable spectator contributions Nf are given by Eq. (51) as 



r + v u 
o 



(74) 



where the effective Wilson coefficients are given by 

°eff,l — °2, <^cfr,2 — °eff,3 — °4j ^ e ff,4 ~ °3, ^eff,5 — °6- 



(75) 



The spectator contributions to the heavy-to-light form factor Ff are given as 



F i = i-^-^ s f B f^m B C^ ti J du^- J dr A 



r+ 



(76) 



where the effective Wilson coefficients are given by 



AT 
C 5 



AT 
C 6 



iV' 



Crfr,4 - C4 + -77, C[ S5 - C 5 + 



(77) 



The final expression can be simplied when we use the definition of /+ given 
in Eq. (59). For example, Tf + Ff is written as 



Tf + Ff = imlCU [(Ci + (1 + F) + f^ 2 F 



4tt AT 

+t—Tra s f 7r f B m B [C 1 + —) J du—— J dr 
im 2 B (U^^C 2 F + im%U{C x + ^) [C(l + K) 



i(r+) 



4vr JV 

Cf /tt/b 



VfUJb f, <f>M f A 
+7ra s / du / dr + 

N m B J u J 



N ■ 



vm 



4tt N 



(78) 



where the definition of /+ in Eq. (59) is used in the last line with mj pa m B = 
2E. And we replace ( by /+ in the term proportional to F. This induces terms 
of 0{a 2 s ), which is neglected. This relation also holds for the combinations 
Tf + Ff for all i. 



In summary, the decay amplitudes af are given by 
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<Z = im%f + fJ(C 1 + %) + ^C 2 F' 



N 



4tt N 



«S = *»i/ + /.[(<* + §) + 



P ■ 2 

a$ = im B , 



a, Cf 



N 



4tt N 



C 3 F' + C 1 (--G(s 1 



p ■ 2 n f j 17 n , Cs\ . ol s C f 
al = im B f + f^l(C 4 + —) + - — 

+ C 3 (^- G(0) - G(l)) + C 4 (^ - 3G(0) - G(s c ) - G(l)) 

+ G 6 (-3G(0) - G(s c ) - G(l)) + G^ 8 (G 5 + C 8 )] |, 

af = -im 2 B f + U[(C 5 + + ^^G 6 (-F' - 12) 



where 



(79) 



(80) 



The decay amplitudes in Eq. (79) at order a s are consistent with the result 
obtained by Beneke et al. [7] at order a s . However, our result goes further 
in the sense that we consider the operators to all orders in a s , as shown in 
Eqs. (64), (68) and (70). At higer orders, our result will be different from the 
result in Ref. [7] since there are two scales m b and /i = \/m b A involved and 
the effects of these two scales cannot be reproduced in the heavy quark limit. 



8 Conclusion 



We have considered the four-quark operators relevant to nonleptonic B decays 
into two light mesons in SCET at leading order in A and to next-to-leading 
order in a s . The construction of the four-quark operators in SCET is process- 
dependent since we first have to specify the directions of the outgoing quarks 
and construct the operators accordingly. In matching onto SCETn, we inte- 
grate out off-shell modes by attaching collinear and soft gluons to fermion 
lines. The result is given as gauge-invariant four-quark operators. The form 
of the gauge-invariant operators is obtained to all orders in a s by using the 
auxiliary field method. The Wilson coefficients of these operators can be com- 
puted by matching the amplitudes between the full theory and SCETi since 
the infrared divergence of the full theory is reproduced in SCETi. And we 
obtain jet functions through the matching between SCETj and SCET n . 
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When the effects of collinear and soft gluons are included, we can obtain gauge- 
invariant operators, and the explicit form of these operators guarantees the 
color transparency at leading order in A but to all orders in a s . Now the idea 
of the naive factorization in which the matrix elements of four-quark operators 
are reduced to products of the matrix elements of two currents has a theoretical 
basis. Furthermore when we include spectator interactions which contribute 
to the nonfactorizable contribution and to the heavy-to-light form factor, the 
amplitudes, which include four-quark and six-quark operators, are factorized 
to all orders in a s . That is, the amplitudes can be written as a convolution of 
short- distance effects represented by the effective Wilson coefficients and long- 
distance effects represented by the light-cone wave functions of mesons. And 
the convolution integrals appearing in the hard spectator interactions and in 
the hard contribution to the form factors are finite. Therefore we have proved 
that the decay amplitudes for nonleptonic B decays into two light mesons at 
leading order in SCET and to all orders in a s are factorized. 



Note that we have not included renormalization group running of the effective 
Wilson coefficients. In order to include the renormalization group running, the 
scaling of the Wilson coefficients can be achieved by calculating the anomalous 
dimensions of the operators, say, from Eqs. (25) and (27). But the running 
of the Wilson coefficients may not be appreciable at this order for the scale 
change from /i = m b to jj, = ^rribA. 



As an application of SCET, we have calculated the decay amplitudes for 
B — > -kit at leading order in SCET. The results are consistent at order a s 
with the result in the heavy quark mass limit, which is shown explicitly here. 
This is not a coincidence because the leading-order decay amplitudes in the 
heavy quark mass limit employing leading-twist meson wave functions corre- 
spond to the leading-order decay amplitudes in SCET. However SCET extends 
the analysis to all orders in a s , and proves that the decay amplitudes are fac- 
torized. The two types of the factorization properties, which correspond to the 
color transparency and the separation of long- and short-distance effects are 
satisfied to all orders in a s in nonleptonic B decays. 



We can go beyond the leading-order calculation and consider subleading cor- 
rections in order to check the validity of the approach using SCET. For exam- 
ple, we can ask questions on how chirally-enhaced contributions can be treated 
in SCET, or how to include higher-twist wave functions of mesons, and how 
we can organize higher-order corrections in SCET. However we stress that this 
is a first step toward understanding nonleptonic B decays, and those questions 
are under investigation. 
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A Derivation of the gauge-invariant four-quark operators using the 
auxiliary field method 

In Section 2 we derived gauge-invariant four-quark operators through the 
matching calculation to order g 2 , in which all the off-shell modes were inte- 
grated out. We can derive the gauge-invariant four-quark operators of the form 
(^Tih) ■ (xr 2 x) using the auxiliary field method. We construct a Lagrangian of 
the on-shell fields with the off-shell fields which are auxiliary fields. Then we 
integrate out the off-shell auxiliary fields to obtain gauge-invariant operators 
entirely in terms of the on-shell fields. In integrating out the off-shell modes, 
the Wilson lines W, W, S and S appear and the final form is gauge invariant. 
The auxiliary field method offers a consistent way to derive gauge-invariant 
operators to all orders. 

The auxiliary field method has been used in SCET for the gauge-invariant 
heavy-to-light current operators [15], and for the processes in which there are 
collinear particles in the n M and n M directions [12]. We extend the idea of the 
auxiliary field method to derive the gauge-invariant four-quark operators in 
nonleptonic decays. Since there are collinear fields both in the n M and the n M 
directions, the auxiliary field method for collinear quarks £ and x is the same 
as the method presented in Ref. [12]. Here we consider how to treat the heavy 
quark interacting with two types of collinear gluons A% and A£. 

For clarity, we first construct the Lagrangian of the on-shell fields with the off- 
shell fields of order p 2 ~ Q 2 , and integrate out the off-shell modes to obtain the 
gauge-invariant four-quark operators. Through this procedure, we obtain the 
Lagrangian and the four-quark operators in SCETi. Once we get familiar with 
the construction of the Lagrangian with the off-shell modes of order p 2 ~ Q 2 , 
we include the off-shell modes of order p 2 ~ QA, and integrate out all the 
off-shell modes. We can obtain the Lagrangian and the four-quark operators 
in SCET n in a gauge-invariant way. 

Let us categorize all the fields in SCETi in powers of A. Note that the small 
expansion parameter A in SCETi is of order \J~aJq and A' in SCETn is of order 
A/Q. In order to avoid confusion, we express the scaling of all the momenta 
in powers of A. The on-shell fields are the collinear fields £, A£, which scale as 
Pn = (n- p,n- p, p±) ~ (A, Q, y/QA), the collinear fields x, A£, which scale as 
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Table A.l 

List of fields to construct the auxiliary Lagrangian in SCETi. The on-shell fields are 
collinear, and heavy quark fields, and the Wilson lines obtained by the corresponding 
gluons are listed in the last column. The off-shell fields, which appear as intermediate 
states, are classified by their momentum scaling behavior. 



Momentum scaling 


Fields 


Wilson lines 


on-shell p% ~ (A, Q, (QA) 1/2 ) 




W 


P £~(Q,A,(QA)V2) 




W 


p£ s ~(A,A,A) 






off-shell K+^~(Q,Q,(QA)V2) 




W Q ,W Q 


Pn ~ (A,Q, (QA) 1 / 2 ) 






^~(Q,A,(QA)V2) 







(Q,A, y/QA), the soft fields q s , A^, which scale as ~ VQA, and the usoft 
fields q us , A^ s , h, which scale as p^ s ~ A. 

In SCETi, we include the following off-shell fields. When the on-shell heavy 
quark h interacts with a collinear gluon A^ (A£), the final heavy quark has mo- 
mentum of order p% ~ (A, Q, y/QR) \pg ~ (Q, A, y/QA)}. We label this off-shell 
heavy quark as ip n ijpn)- If the on-shell heavy quark h interacts with A^ and A£, 
the final off-shell heavy quark has momentum of order p^+Pjt ~ (Q, Q, VQA), 
and p 2 ~ Q 2 . We label this off-shell field as ipn- A gluon interacting with A^ 
and A£ through the triple gluon vertex also has a momentum + p£, and we 
label this gluon as Aq. We list all the on-shell fields and the off-shell auxil- 
iary fields in Table A.l. The collinear quarks can also be off-the mass shell of 
order p 2 ~ Q 2 , which we denote as £q and \Q- The treatment of the off-shell 
collinear quarks is presented in Ref. [12], and we will not repeat it here. Since 
the soft and usoft quarks and gluons are not relevant in our argument, we do 
not present them in Table A.l. 

The auxiliary gluon Lagrangian in SCETi is given by 
where the covariant derivative Dq is given by 



2 



■n 



iD% = —(n-V + gn-A n ) + —(n-Q + gn- A*) 

+V1 + gA^ + Ql + gA^, (A.2) 

and only the leading terms in Eq. (A.2) is included in £| ux at leading order. 
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We can separately obtain the solution of the first term and the second term 
in Eq. (A.l). The equation of motion of the first term is given by 



[iD Qli + gA Q ^ [iDq + gA Q , iD Q + g A%\] = 0. 

The leading-order solution is obtained by making an ansatz 
WqWq = WW\ 



(A.3) 



(A.4) 



Here Wq and Wq are essentially the Fourier transform of the Wilson lines 





r y 




W Q (y) =Pexp< 


{"J 


ds 




*■ — oo 






r y 




W Q (y)=Pexp< 


* / 


ds 




*■ — oo 





which satisfy 

(V + gA Q + gA n )\V Q = 0, n ■ (Q + gA Q + gA^Wq = 0. 



n 



(A.5) 



(A.6) 



The Lagrangian with auxiliary heavy fields is complicated because only the 
sum of three graphs shown in Eq. (16) is simple, but the individual diagrams 
show complex behavior. Recall that Eq. (16) is given by 

n/r i n/r ^9 2 f -T-, rp n ' A^n ■ A c 2 n ■ A n n ■ A Tl 
M a + M b = —fabcq^iTa - - - b + g qT 1 - - — b, (A.7) 



n-qjtn-q n 



n-q n n-qn 



and we construct the Lagrangian such that the first term in Eq. (A.7) is 
reproduced by the auxiliary field attached with the triple gluon vertex with 
A^ and A£, and the second term is produced by attaching A^ and A^ starting 
from the vertex in this order. It means that when a heavy quark h becomes an 
off-shell field ipn interacting with A 1 ^ and the heavy field h first interacts 
with A£ to become ife, and then interacts with A^ to become ipn- We can 
choose a different convention in which h interacts with A^ first and then 
interacts with A£ to become ipH- Both methods give the same result, as will 
be shown later. However, it is necessary to specify the order of the interaction 
since the Lagrangian for the off-shell fields is obtained by expanding the QCD 
Feynman rules in powers of A. 

With the prescription described above, the Lagrangian for the heavy quark 
with the auxiliary fields is given as 
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£lx = *l> n 9n ■ A nh + • A nh + ip H gn ■ A Q (h + lf) n ) 
+i) H gn ■ (A Q + A n )ipn + h.c. 
+^ n (n -V + gn- A n )ijj n + ?/%(n -Q + gn- A„)^n 



+i) H [n ■ V + gn ■ (A Q + A n )] ^ H 



(A.8) 



Note that there is no term in which ip n interacts with A£ to become ipn 
according to our prescription. The inclusion of the auxiliary field Aq in the 
second line of Eq. (A.8) should be verified at order g 3 , but these terms are 
included since they are kinematically allowed. 

Solving for ipn, ipn and ipn from Eq. (A.8), we obtain 



gn- A^h + n- (Q + gAn)ipn = 0, (A.9) 
gn ■ A n h + n ■ (V + gA n )^ n = 0, 

gn ■ A Q (h + i/j n ) + gn ■ (A n + A Q )fa + n ■ \v + g(A n + A Q )U H = 0. 



The first equation in Eq. (A.9) can be solved for ■?/% as 

fa=(W-i)h, 



(A.10) 



and adding the second and the third equations in Eq. (A.9) yields 

if> n + i(>H = (W Q -l)(h + fa). (A.ll) 



Therefore the heavy quark field can be written as 

h + V% + V>n + 1>H = h + fa + (Wq - l)(h + fe) = W Q Wh. 



(A.12) 



We can arrive at an equivalent conclusion by rewriting Eq. (A. 7) as 

iQ 2 r - ^ n • A%-n - At. , _ n ■ A Tl n ■ A n 

M a + M b = ^-fa^qT^a 2= + g 2 qT 1 -=^, 

2 n-qrjn-q n n ■ q^ nb ■ q n 



(A.13) 



where the order of the gluons in the second term is reversed. In this case, we 
require that the heavy quark h interacts with A^ first and then interacts with 
A£ to generate ipn- The corresponding Lagrangian is given by 



£lx = ^n9n ■ A n h + ip w gn ■ A w h + ip H gn ■ A Q (h + fa) 
+tp H gn ■ (A Q + An)ip n + h.c. 
+^ n (n -V + gn- A n )ip n + ip w (n -Q + gn- A^i/% 
+ip H n- Q + gn - (A Q + A w ) tp H - 



(A. 14) 
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Solving the equations of motion, we obtain 



h + 1> n + fa + iPH = WQWh = W Q Wh, (A.15) 



where the last equality comes from the ansatz WqWq = WW . This is equiv- 
alent to the result in Eq. (A. 12). 

Now we construct the Lagrangian including the off-shell modes of order p 2 ~ 
QA. Through this procedure, we obtain the Lagrangian and the four-quark 
operators in SCETn. In SCETn, we have to integrate out all the off-shell modes 
in Table A. 2 by constructing the Lagrangian with additional auxiliary fields 
with additional soft momentum, with the index X. The on-shell and the off- 
shell fields in SCET n are listed in Table A. 2. Note that the momentum scaling 
of each fields has changed in order to accomodate SCETn- Since Pn+P^+Pg ~ 
p^ + pft, we do not put the index X to ipn and Aq. For the same reason, we 
express the off-shell heavy quarks as ip n and 

We include the off-shell modes due to the soft momentum. We label these 
fields with the index X to denote the off-shellness p 2 ~ QA. The auxiliary 
Lagrangian including the off-shell modes of order p 2 ~ QA is given by 



£aux = i>n9~n ■ ( A n + A n)h + ^9 n ' ( A W + A n)h + ^ H gfl ■ A Q (h + 

+^ H gn ■ (A Q + A* + A n )^ + h.c. 



+^ n n ■ [V + g(A$ + A n )\ ^ n + ■ \ Q + g(A£ + A„)\ ^ 
+$ H n ■ [V + g(A Q + A n )] fa. (A. 16) 

Note that there is no soft gluon A^ involved in Eq. (A. 16) for the heavy quark. 
Table A.2 

List of fields to construct the auxiliary Lagrangian in SCETn. The fields with 
the index X respresent the off-shell fields with the soft momentum. The momen- 
tum scaling of ipn, Aq, ip n and ipn does not change with the addition of the soft 
momentum. 



Momentum scaling 


Fields 


Wilson lines 


on-shell p% ~ (A 2 /Q,Q,A) 


e, a% 


W 


p^(Q,A 2 /Q,A) 


X,A» 


W 


p£~(A,A,A) 


q s , A%, h 


S, S 


off-shell p£+p£ + Ps~(Q,Q,A) 


*Ph, A» Q 




p£ + p^~(A,Q,A) 




W x , S x 


< + l^~(Q,A,A) 


^n, XX, Ait x 


Wx,S x 
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It appears in the Lagrangian for the collinear quarks. The explicit Lagrangian 
for the collinear quarks involving the off-shell fields of order p 2 ~ QA can be 
found in Eqs. (A4) and (A5) in Ref. [12]. 

The difference between the Lagrangians given in Eqs. (A. 8) and (A. 16) is that 
all the off-shell fields are labeled by X, and A n — > A* + A n , A„ — > A^ + A„. 
Therefore the heavy quark field, after the off-shell fields are integrated out, is 
written as 



h + i> n + Ipn + i>H = W$W X h, 



(A.17) 



and the collinear quarks are written as [12] 



-x 



£ + + & = W Q S X Z, x + Xx + Xq = W^S x x, 



(A.18) 



where Wq , and Sx are the Fourier transforms of the Wilson lines 



W Q (y) = P exp< ig / ds n ■ A Q (sn) + n ■ A nX {sn) + n ■ A n (sn) 



oo 

y 



S x (y) = Pexp 



"J 



ds 



n ■ A nX {sn) + n ■ A s (sn) 



(A.19) 



and Wq and S x are obtained by replacing by n M in Eq. (A.19). The Wilson 
lines satisfy 



W^Wg = W X W X , S X W X = WS\ S X W X = WS ] . 



(A.20) 



The last two relations in Eq. (A.17) were obtained in Ref. [12] and the first 
relation is new. 

By integrating out the off-shell fields, the singlet four-quark operator in SCET n 
is given by 



X (x + Xx + Xq) T 2(x + Xx + Xq) 
IS^W^YAW^Wxh) ■ (xS x W^F 2 W^S xX 



= (^Wx^w^Wxh) ■ ( X T2x) = twrt&h ■ xr 2 x, (A.21) 

using Eq. (A.20), and similarly the nonsinglet four-quark operator is given as 
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= (tsWx)r 1 (s^W x h) a -x a T2(W j x s x x) p 

= (zws^r^h)* ■ {xW) a V 2 {W\)p. (A.22) 

This result is the same as the explicit calculation obtained in Section 2, and 
it is a proof to all orders in a s using the auxiliary field method. 



B Derivation of the subleading operators using the auxiliary field 
method 

We can also derive the subleading operators ( j la ' lb ^ in Eq. (36) and J- la ^ and 
4 16) in Eq. (53) using the auxiliary field method. Here we have to consider 
off-shell modes from the collinear gluons both in the n M and n M directions. For 
simplicity, we consider the subleading operators in SCETj disregarding the off- 
shell modes by soft gluons. The leading-order result was derived in Ref. [12], 
and here we present a new result which yields gauge-invariant operators at 
subleading order. In addition to the solutions in Eq. (A. 3) at leading order, 
when we include the subleading terms in the covariant derivative in Eq. (A. 2), 
we obtain the equations of motion at subleading order, which are given as 

[W Q n ■ VWl [W Q n ■ QWl iV» Q± }\ = 0, 

Wqn ■ QWl [W Q n ■ VW^ iP$J] = 0, (B.l) 
where iT>j^ ± = iD^ L + gA^ ± . Here we make an ansatz 

W^V: ±Q = iD» nA Wl or W^ ±Q = iD^Wl (B.2) 
and they satisfy Eq. (B.l). 

First the intermediate form of the subleading four-quark operators relevant to 
the nonfactorizable spectator contribution and the heavy-to-light form factor 
can be obtained by neglecting the off-shell modes in the n M direction since we 
are interested in the operators proportional to gA£ at leading order in g. As 
a definite example, let us concentrate on the operator and obtain the 
gauge-invariant form by integrating out the off-shell modes. The derivation 
for the operator can be done in a similar way. In SCETi, without the 

off-shell modes in the n M direction, the operator of^ is given by 
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x (x + x )r 2i K^H) + ^Qj|(x + x Q ), (b.3) 

where the color indices are suppressed. When we include off-shell modes, each 
field is written as 

h + ^ n + *Pn = W Q Wh, £ + i Q = W Q ^ X + XQ = W QX . (B.4) 



Then the operator in Eq. (B.3) is written as 



{iw ] Q )Y u {w Q wh)\ ■ [(xW^r 2i w Q ^w* Q ([ip n±Q w Q }^x) 



(B.5) 



If we use the ansatz in Eqs. (A. 4) and (B.2), the singlet operator becomes 



u s 



^wUw Q ) a Ti(Wh) a 



x 



mh )/»r 2 {Wq^—qWq lip n±Q w Q }^x) 



2 / V/3. 



(W Q rA 



X 



{xW ] Q w Q ) p v 2 — -([ijp^wlWQY-x) 



2^/3. 

1 ^„ + .^ t 



(B.6) 



and the nonsinglet operator becomes 



og b) = [(^)^r 1 (w w/i) Q ] • [(xwl) a r 2 (w Q ^wl[ip n±Q w Q }^x) 



'2'V/3 
'2^/3 



;bj) 



These are the sub leading operators which are gauge invariant. Other operators 
can be shown to have the forms presented in the paper in a similar way. 
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